
lapop12 = readRDS('data/lapop12.rds')

lapop12 = lapop12 %>% 
  mutate(
    soc.cl_missing = as.numeric(is.na(soc.cl_bin) == TRUE, 1,0), 
    torture_bin_missing = as.numeric(is.na(torture_bin) == TRUE, 1,0), 
    fair_trial_missing = as.numeric(is.na(fair_trial) == TRUE, 1,0)
  )


m1 = lm_robust(soc.cl_missing ~ trial + q1 + q2, 
               data = lapop12,
               fixed_effects = estratopri + estratosec)

m2 = lm_robust(torture_bin_missing ~ trial + q1 + q2, 
               data = lapop12,
               fixed_effects = estratopri + estratosec)

m3 = lm_robust(fair_trial_missing ~ trial + q1 + q2, 
               data = lapop12,
               fixed_effects = estratopri + estratosec)
texreg::texreg(list(m1, m2, m3), include.ci = F, 
               stars = c(.001, .01, .05, .1), 
               symbol = '\\dagger')
